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Abstract:- In this work, forced convection heat transfer inside a channel of rectangular 
section, containing some rectangular baffle plates, is numerically analyzed. We have 
developed a numerical model based on a finite-volume method, and we have solved the 
coupling pressure-velocity by the SIMPLE algorithm [6]. We show the effects of various 
parameters of the baffles, such as, baffle's height, location and number on the isotherms, 
streamlines, temperatures distributions and local Nusselt number values. It is concluded 
that: I) the baffles location and height has a meaningful effect on isotherms, streamlines 
and total heat transfer through the channel. (II) The heat transfer enhances with 
increasing both baffle's height and number. 
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I. INTRODUCTION 

In recent years, a large number of experimental and numerical works were 
performed on turbulent forced convection in heat exchangers with different type of 
baffles [1-4]. This interest is due to the various industrial applications of this type of 
configuration such as cooling of nuclear power plants and aircraft engine ... etc. 
S.V Patankar and EM Sparrow [1] have applied a numerical solution procedure in order to 
treat the problem of fluid flow and heat transfer in fully developed heat exchangers. 
These one was equipped by isothermal plate placed transversely to the direction of flow. 
They found that the flow field is characterized by strong recirculation zones caused by 
solid plates. They concluded that the Nusselt number depends strongly on the Reynolds 
number, and it is higher in the case of fully developed then that of laminar flow regime. 

Demartini et al [2] conducted numerical and experimental studies of turbulent flow 
inside a rectangular channel containing two rectangular baffles. The numerical results 
were in good agreement with those obtained by experiment. They noted that the baffles 
play an important role in the dynamic exchangers studied. Indeed, regions of high 
pressure 'recirculation regions' are formed nearly to chicanes. 

Recently, Nasiruddin and Siddiqui [3] studied numerically effects of baffles on forced 
convection flow in a heat exchanger. The effects of size and inclination angle of 
baffles were detailed. They considered three different arrangements of baffles. They 
found that increasing the size of the vertical baffle substantially improves the Nusselt 
number. However, the pressure loss is also important. For the case of inclined baffles, 
they found that the Nusselt number is maximum for angles of inclination directed 
downstream of the baffle, with a minimum of pressure loss. 

More recently, Saim et al [4] presented a numerical study of the dynamic behavior of 
turbulent air flow in horizontal channel with transverse baffles. They adapted numerical 
finite volume method based on the SIMPLE algorithm and chose K-s model for 

treatment of turbulence. Results obtained for a case of such type, at low Reynolds 
number, were presented in terms of velocity and temperature fields. They found the 
existence of relatively strong recirculation zones near the baffles. The eddy zones are 
responsible of local variations in the Nusselts numbers along the baffles and walls. 

We know that the primary heat exchanger goal is to efficiently transfer heat from 
one fluid to another separated, in most practical cases, by solid wall. To increase heat 
transfer, several approaches have been proposed. We can cite the specific treatment of 
solid separation surface (roughness, tube winding, vibration, etc.). This transfer can also 
be improved by creation of longitudinal vortices in the channel. These eddies are 
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produced by introducing one or more transverse barriers (baffle plates) inside the 
channel. The formation of these vortices downstream of baffles causes recirculation zones 
capable of rapid and efficient heat transfer between solid walls and fluid flow. It is this 
approach that we will follow in this study. Indeed, we are interested in this work on the 
numerical modeling of dynamic and thermal behavior of turbulent forced convection in 
horizontal channel where two walls are raised to a high temperature. This channel may 
contain one or several rectangular baffles. 

A special interest is given to the influence of different parameters, such as height, 
number and baffle positions on heat transfer and fluid flow. 

II. MATHEMATICAL FORMULATION 

The geometry of the problem is shown schematically in Figure 1. It is a rectangular 
duct with isothermal horizontal walls, crossed by a stationary turbulent flow. The 
physical properties are considered to be constants. 
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Figure 1: The studied channel 



At each point of the flow the velocity has components (u, v) in the x and y 
directions, and the temperature is denoted T. The turbulence modeling is handled by the 
classical model (k-s). k is the turbulent kinetic energy and s the viscous dissipation of 
turbulence. 

The Transport equations (continuity, momentum, temperature, turbulent kinetic 
energy and dissipation of turbulence) governing the system, are written in the following 
general form : 
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Where p is the density of the fluid passing through the channel and and are given 

by: 
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The constants used in the turbulence model (k-£) Are those adopted by Chieng and 
Launder (1980) [8] : 
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Boundary conditions 

At the channel inlet: 

u = U in ;v = and T = T in 

k = 0,005 Uf n and s = 0,1 k 3/2 
At solid walls: 

u = v = \k — s — and T = T w >- T in 

At the channel exits the gradient of any quantity, with respect to the longitudinal 
direction x is nul. 

Our goal is to determine velocity and temperature fields, as well as the turbulence 
parameters. Particular attention is given to the quantification parameters reflecting heat 
exchanges such as the Nusselt number (local and average). 

III. NUMERICAL FORMULATION 

The computer code that we have developed is based on the finite volume 
method. Computational domain is divided into a number of stitches. To choose the number 
of cells used in this study, we performed several simulations on a channel (with or 
without baffles). Finally, we have opted for a 210 x 90 meshes. 

The mesh dimensions are variables; they tighten at the solid walls neighborhoods. 
Consequently, the stitch density is higher near the hot walls and baffles. 

We consider a mesh having dimensions Ax and Ay. In the middle of each volume 
we consider the points P, called centers of control volumes. E, W, N, S are the centers of 
the adjacent control volumes. We also consider centers, EE; WW, NN, SS. The faces of 
each control volume are denoted e, w, n, s. 
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Figure 2: Mesh with P at center 



By integrating the transport equation (1) of $ on the control volume, 
relation to the x direction: 
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Where 



<j j, is the source term 

L <P Ip 



I w 



One note F = pu , the convection flow and D = T/Ax, the diffusion coefficient. And 
with taking : 

F e = [pu] e \F w =\pu\ 

D e =[r/Ax] e ;D w =[T/Ax\ 

And : [d(f\ e =<f> E -0 P ; [d(/)] w =0 P ~0 W 
We found : 

p e & -F w <t> w = D e (</> E - - d w {$ p - 4) + [sX 

To estimate $ on the faces " e " and " w " we opted for the classical quick scheme 
[5] which is a quadratic forms using three nodes. The choice of these nodes is dictated by 

the direction of flow on these faces (u >- or u -< 0) . 

Finally, the transport equation (1) is discretized on the mesh with P at center as : 

*P <f>P = OW K + a E <I>E + a WW Kw + a EE <I>EE + [ S \ P ( 2 ) 

With : 

6 13 
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The same thing is done for the vertical "y" direction, using the faces "n" and "s" 
and by introducing the Quick diagram nodes N, NN, S and SS. 

if convection flow and diffusion coefficient of the sides "e", "w", "s", "n" are known 
And especially the source term S ± , the solution of equation (2) then gives us the values 

of $ in different P nodes. 

One notes that to accelerate the convergence of equation (2) we have introduced a 

relaxation factor : 

a 



p 



(j) p = % + a E (j) E + % W <f>WW + a EE <f>EE + [S A P + a P<f>P ( 3 ) 

Where (j)p is the value of (j) p in the previous step. 

The source term appears especially in the conservation of momentum equations in 
the form of a pressure gradient which is in principle unknown. To get around this 
coupling we have chosen to use in our code the "Simple" algorithm developed by 
Pantankar [7]. The basic idea of this algorithm is to assume a field of initial pressure and 
inject it into the equations of conservation of momentum. Then we solve the system to find 
a field of intermediate speed (which is not fair because the pressure isn't.) The continuity 
equation is transformed into a pressure correction equation. This last is determined to find 
a pressure correction that will inject a new pressure in the equations of motion. The 
cycle is repeated as many times as necessary until a pressure correction equal "zero" 
corresponding to the algorithm convergence. In the end we solve the transport equations 
of T, k and £. 

In this approach, a problem is encountered. It is known as the checkerboard 
problem. 

The risk is that a pressure field can be highly disturbed by the sensed formulation 
which comprises performing a linear interpolation for estimate the pressure value on the 
facets of the control volume. To circumvent this problem we use the so-called staggered 
grids proposed by Harlow and Welch [6]. In this technique, a first grid pressure (and other 
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scalar quantities T, k and £) is placed in the center of the control volume. While other 
staggered grids are adopted for the velocity components u and v (see Figure 3). 
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Figure 3: Shifted mesh 



The scalar variables, including pressure, are stored at the nodes (I, J). Each node 
(I, J) is surrounded by nodes E, W, S and N. The horizontal velocity component u is stored 
on faces "e" and "w", while the vertical component v is stored on the faces denoted "n" 
and "s". 

So the control volume for pressure and other scalar quantities T, k and £is 

[(/ + 1, j); (i 1); 0\ j + 1); (U j)] 

For component u control volume is 

[(/,;);(/, j + -1,7 + 1);(/ 

While for the v component we use: [(l+l, / -1); (i + 1, /);(/,/);(/ — 1, / — 1)]. 
By integrating the conservation momentum equation in the horizontal direction on the 

volume [(/, 7); (7,7 + 1); (7-1, 7 + 1); (7-1, 7)]. w e found : 

<h,Pu =Z«A + [^-U)-P(/,/)h +1 -y.) (4) 

Where : ^^ci nb u nb = a w w w + a w u w + a EE u EE + a E u E + a S s u ss + a s u s + a NN u NN + a N u N 

The coefficients a^are determined by the quick scheme mentioned above. 

Similarly, the integration of conservation momentum equation in the vertical 
direction on the volume 

[(/ + 1, J - 1); (i + 1, J); (i, J); (i - 1, J - 1)], 
gives us : a, , jVl ,j =£a nb v nb +[P(I,J-l)-P(I,J)](x i+1 -Xj) (5) 

One considers primarily an initial pressure field P . The provisional solution of the 
equations (4) and (5) will be denoted u and v . We note that u and v doesn't checks the 
continuity equation, and: 

a u u *i.J = ^2l a nb u *"b +[P\l-\,J)-P\l,J)\(y j+l -yj) (6 -a) 
a IJ v\j=^aJ'rt+[p\l,J-l)-P\l,J)](x i < 1 -x t ) (6-b) 

At this stage any one of the three variables is correct. They require correction : 

'p'=P-P* 

u - u - u (7) 



V = v - V 



Injecting (7) in equations (4), (5) and (6) we find : 

a u u' u = ^ nb u\ b +[F (7 - 1, J) - P (I, J)](y j+l - y, ) (8 - a) 

V'u = +[ p ' ( 7 > J - D - F ( 7 ' - x i ) ( 8 - b) 
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At this level, an approximation is introduced. In order to linearize the 
equations (8), the terms ^ \a nh U \ h an d^£^v'^ are simply neglected. 

Normally these terms must cancel at the procedure convergence. That is to say that 
this omission does not affect the final result. However, the convergence rate is changed 
by this simplification. It turns out that the correction P' is overestimated by the Simple 
algorithm and the calculation tends to diverge. The remedy to stabilize the calculations is 
to use a relaxation factor. 

We Note that a further treatment of these terms is proposed in the so-called 
algorithms "SIMPLER" and "SIMPLEC" [7]. 
The equations (8) become: 

u\j = d u [F(I -I, J)- F(I, J)] (9 - a) 
v'u = d^FihJ-V-FihJ)} (9-b) 

with du = y J +1 ~ y i and di . = ~ x > 

Equations (9) give the corrections to apply on velocities through the formulas (7). We 
have therefore : 

u u = uu + d u [F (I-1,J)-F (I, J)] (10 - a) 
v Lj = vV,+djF(/,/-l) -/*(/,/)] (10-6) 

Now the discretized continuity equation on the control volume of scalar quantities is 
written as follows : 

[{puA) M J - (puA) u J- [(pvA) I j+l - (pvA), . J = 

Where an are the size of the corresponding faces. 

The introduction of the velocity correction equations (10) in the continuity 
equation gives a final equation allowing us to determine the scope of pressure corrections 
P: 

a\jF(I,J) = a\ +lJ P' (/ + W) + a\_ u F (I - 1, /) + a\ J+x F (/, / + 1) + d n _ x F (/, J -I) (11) 

The algorithm can be summarized as follows : one starts from an initial 
field P*,u*,v* and <j) , with (j) represents the scalar quantities T, k and 8. Then the system (6) 

is solved to have new values of u and v . 

Then the system (11) is solved for the corrections field pressure P'. 

Thereafter, pressure and velocity are corrected by equations (7) and (10) to have P, u and 
v. 

We solve the following transport equation of scalar quantities (2), ^ = T, K and £ . 
Finally, we consider: 

P = P ;u —u \v* =v\(/)* 

And the cycle is repeated until convergence. 

IV. RESULTS AND DISCUSSIONS 

The dimensions of the channel presented in this work are based on experimental 
data published by Demartini et al [2]. The air flow is carried out under the following 
conditions. 

> Channel length: L = 0.554 m; 

> Channel diameter D = 0.146 m; 

> baffle height : <h <0 1. m; 

> baffle thickness: 8 = 0.01 m ; 

> Reynolds number: Re = 8.73 10 4 ; 

> The hydrodynamic and thermal boundary conditions are given by 
At the channel inlet 

u = u in =7,8m/s; 



v =o ; r = 7;=3oo K 
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On the channel walls 

a = v = and T = T w = 313 K 

At the channel exit the system is assumed to be fully developed, ie 




O 



First, we compared the structure of streamlines and isotherms in a channel without 
baffle with those obtained when we introduce baffle having a height h = 0.05 

m (Case 2-1) at the abscise L x = 0.218 m. Indeed, in figures 4 and 5 we present streamlines 
and isotherms for two configurations. 

These results clearly show the importance of presence of baffle (acting as a 
cooling fin). Indeed the baffle increases heat transfer between the wall and the fluid. This 
increase is caused by recirculation zones downstream of the baffle. 



(a) 





1 







(b) 

Figure 4: streamlines: (a) case 1; (b) case 2-1 
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(a) 

















(b) 

Figure 5: Isotherms: (a) case 1; (b) case 2-1 



We then studied the influence of baffle height on flow structure and heat 
transfer. We chose the heights (h = 0.05 m; cases 2-1) h = 0, 073 m (Case 2-2) and h 
= 0, 1 m (case 2-3). 

The figures 6 and 7 present streamlines and isotherms for a channel containing baffle with 
different heights. 
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It is clear that the isotherms are more condensed near baffle, as and when the 
height h increases. This indicates an increase of heat transfer near baffle. Indeed, the 
increase in h expanded exchange area between fluid and walls. In addition, when h 
increases, the recirculation zone becomes increasingly important (see Figures 4 (b) and 
figures 6 (a) and (b)). This causes an acceleration of flow, which improves heat transfer 
within the channel. 

We present figures 8 and 9 in order to identify the influence of h on velocities and 
temperature profiles along y at x = 0.45 from channel inlet. Four values of h are 
considered. 
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31 6 aid- 
Figure 9: Temperature profiles along y at x = 0.45 m « influence of h 



» 



These results allow us to conclude that the increase of baffle height has two 
contradictory effects. It is true that there is substantial increase in heat exchange 
(appearance of recirculation zones most common), but there is also a loss of pressure 
(flow blockage). 

In the purpose of measuring the influence of 'h' on local heat transfer, we have 
presented in figure 10, the local Nusselt number along the channel for the four cases of h 
considered. It shows that upstream of the baffle (x <0.2 m) curves are confused , while, 
just after baffle location, effect of baffle height on Nusselt number become increasingly 
important away from the baffles. 
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Figure 10 : Local number Nu along the channel « influence of h » 

We then studied the influence of baffles number and position 'Figure 11 and 12'. 











(a) 









(b) 
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(c) 

Figure 11 : Streamlines : (a) cas 3-1; (b) cas 3-2; (c) cas 3-3 

We consider two baffles of same height h = 0.08 m. The first baffle is placed at the 
distance L x = 0.15 ra, while the second is positioned at the distance d i = 0.05 m from the 
first (case 3 -1), di = 0.10 m (case 3 -2) and di = 0.15 m (case 3 -3). We note the 
existence of two recirculation zones downstream of the first baffle. Also the first 
recirculation zone 'defined by baffles' becomes increasingly important, which contribute 
to an increase in heat transfer in this area, as shown in Figure 12. Indeed, when the 
baffles are close to each other 'di small', the fluid is blocked in the chimney delimited by 
x = Li and x = Li + di. This reduces the flow speed at that location and thus there will be 
a decrease in the heat transfer in this area. Or when di increases, the fluid has sufficient 
space to move rapidly, hence heat transfer increases in this zone 'see figurel5 '. 




(c) 

Figure 12 : Isotherms (a) : « cas 3-1 » ; (b) : « cas 3-2 » ; (c) « cas 3-3 » 

In the following figures 13 and 14, we compare the profile of horizontal velocity 
and temperature distribution. Calculation is done on the y-axis at x = 0.45 m from the 
channel inlet for three different baffles spacing. We find that more spacing is important 
more heat exchange is important in areas limited by two baffles. 
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Figure 13 :Profiles of the horizontal velocity at x = 0.45 m « Influence of the spacing 

between two baffles » 
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Figure 15 present the local Nusselt number along the channel for three considered 
spacing d i = 0.05m, 0.1m and 0.15m. 
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Figure 15 :Local Nusselt number along the channel « Influence of baffles number and 

position » 
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We distingue two areas : 

The first defined by <x <0.3m while the second by x> 0.3m. For the first region, 
an increase in number of baffles improves the heat transfer; whereas the reverse is true 
for the second zone. 

We chooses to show the influence of baffles number on heat transfer along the 
channel for three cases: Channel without baffles, channel containing a single baffle and 
channel containing two baffles 'Figure 16'. 

It should be noted that generally, heat transfer is proportional to baffles number. 
Indeed, the Nusselt number characterizing heat transfer within the channel increases with 
increasing baffles number. 



_ c c 




_ 0,2 0,4 0,6 

Figure 16: Local Nusselt number along the channel 

V. CONCLUSION 

The thermal behavior of a stationary turbulent forced convection flow within a 
baffled channel was analyzed. The results show the ability of our code to predict dynamic 
and thermal fields in various geometric situations. We studied mainly the influence of 
baffles height and spacing on heat transfer and fluid flow. One can conclude that: 

1 - Increase in the baffle height improves heat transfer between channel walls and 
fluid passing through it; 

2 - Heat transfer becomes increasingly important with adding baffles. 

3 - The spacing di between baffles has different effects on local heat transfer. Any 
time di has not a lot of influence on the overall heat transfer in the channel. 

In perspective, we intend deepen and clarify our results. Indeed, we will adapt our code to 
others geometric cases (non-rectangular baffles or baffles inclined). Finally, we will also 
try to refine more the turbulence model. 
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